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Abstract 



We consider here a generalization of the random quantum rotor model in 
which each rotor is characterized by an M-component vector spin. We focus 
entirely on the case not considered previously, namely when the distribu- 
tion of exchange interactions has non-zero mean. Inclusion of non-zero mean 
permits ferromagnetic and superconducting phases for M = 1 and M = 2, 
respectively. We find that quite generally, the Landau theory for this system 
can be recast as a zero-mean problem in the presence of a magnetic field. Nat- 
urally then, we find that a Gabay- Toulouse line exists for M > 1 when the 
distribution of exchange interactions has non-zero mean. The solution to the 
saddle point equations is presented in the vicinity of the bi-critical point char- 
acterized by the intersection of the ferromagnetic (M=l) or superconducting 
(M=2) phase with the paramagnetic and spin glass phases. All transitions 
are observed to be second order. At zero temperature, we find that the fer- 
romagnetic order parameter is non-analytic in the parameter that controls 
the paramagnet/ferromagnet transition in the absence of disorder. Also for 
M=l, we find that replica symmetry breaking is present but vanishes at low 
temperatures. In addition, at finite temperature, we find that the qualitative 
features of the phase diagram, for M = 1, are identical to what is observed 
experimentally in the random magnetic alloy LiHo x Y\- x F\. 

I. INTRODUCTION 

Transport in granular metals is mediated by activated transport among the metallic 
grains. In granular superconductors composed of spatially separated metallic grains, such 
single particle charging events ultimately lead to a destruction of phase locking between the 
grains. This state of affairs obtains because the particle number, n, and phase, 9, associated 
with each grain are conjugate variables. If two grains phase lock, the resultant infinite 
uncertainty in the particle number leads necessarily to single particle charging. Should the 
single particle charging energy sufficiently exceed the Josephson coupling energy between 
grains, superconductivity is quenched^. 

The simplicity of the physics underlying the quantum phase transition in an array of 
superconducting islands implies that the resultant Hamiltonian, 
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H=~EcE{^) -E^cos^-^), (1) 

is characterized by only two parameters: 1) a charging energy, Ec, and 2) a Josephson 
coupling energy, J it j. In writing Eq. ([]]), we assumed that the islands occupy regular sites 
on a 2D lattice and only nearest neighbor, Josephson coupling is relevant. For the 

ordered case in which = J , the superconductor- insulator transition is well-studiedi as 
the parameter E c /Jq increases. In sufficiently disordered systems, however, the Josephson 
energies are not all equal and in fact can be taken to be random. 

We are concerned in this paper with the case in which the Josephson energies are random 
and characterized by a Gaussian distribution 



p(J " )= 7ib exp 



(Jjj - Jo)' 
2 J 2 



(2) 



with non-zero mean, J . While random Josephson systems have been treated previouslyEHl, 
such studies have focused predominantly on the zero mean case in which J = 0. This 
limit differs fundamentally from the non-zero mean case, because the superconducting phase 
exists only when J 7^ 0. Hence, for J — 0, there is an absence of an ordering transition. 
Nonetheless, the zero-mean case is still of physical interest because as Read, Sachdev, and 
Ye (RSY)i as well as Miller and HuseS have shown, a zero temperature transition from a 
quantum spin-glass to a paramagnet occurs as the strength of the quantum fluctuations 
increases. To put the random Josephson model in the context of the work of RSY, it is 
expedient to introduce the change of variables 

Sj = (cos 9i, siia 9 i). (3) 

The resultant Josephson Hamiltonian 



H =- E oE[»a) (4) 




is recast as a 2-component (M=2) interacting spin problem with random magnetic interac- 
tions. This model is easily generalizable to describe interactions among any M-component 
spin operator (or quantum rotor) in the group O(M). M=l corresponds to Ising spins and 
is relevant to random magneticQ systems such as LiHoo.i67Yo.833F4 whereas the M=3 limit 
is applicable to spin fluctuations in the cupratesi. In this paper we focus primarily on 
the M=l and M=2 cases in which the ordered phases correspond to a ferromagnet and a 
superconductor, respectively. 

Two types of fluctuations control the transitions between the phases: 1) dynamic quan- 
tum fluctuations arising from the charging energy and 2) static fluctuations induced by the 
disorder. For non-zero mean, three phases (spin-glass, paramagnet, ferromagnet (M=l) or 
superconductor (M=2)) are expected to meet at a bi-critical point. Experimentally, a bi- 
critical point is anticipated whenever J , J and Ec are on the same order of magnitude. It is 
the physics at this bi-critical point that we focus on in this paper. To discriminate between 
these phases, we distinguish between thermal averages, (...) and averages over disorder, [...]. 
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In the superconductor (M=2) or ferromagnetic phases (M=l), the disorder and thermal av- 
erage of the local spin operator, [(«%„)] 7^ 0, is non-zero. In the spin glass phase, the thermal 
average (S iu ) 7^ 0, while [(«Si„)] = 0. At zero temperature, quantum fluctuations and static 
disorder conspire to lead to a vanishing of the static moment in the paramagnetic phase, 
that is, (Si V ) = 0. RSY have performed an extensive study of the spin-glass/paramagnetic 
boundary using the replica formalismi in the case of zero mean in which the purely ordered 
phase is absent. We adopt this formalism here in our analysis of the bi-critical region. The 
only prior study on the non-zero mean case is that of Hartman and WeichmanEl who studied 
numerically the spherical limit, M — ► 00, of the quantum rotor Hamiltonian. In their study, 
they found that the spin- glass phase is absent in the M — > 00 limit for d = 2. In the present 
work, we will not address the issue of dimensionality because we limit ourselves to a mean- 
field description in which all fields are homogeneous in space. Our results then are valid 
above some upper critical dimension that can only be determined by including fluctuations 
around the mean-field solution. 

This paper is organized as follows. In section II, we use the replica formalism to average 
over the disorder explicitly and obtain the effective Landau free-energy functional. We show 
that the leading terms in the Landau action resulting from the non-zero mean are analogous 
to those arising from an external magnetic field in the zero- mean problem. As a result, the 
appropriate saddle-point equations can be solved using a direct analogy to the zero-mean 
problem in the presence of a magnetic field. Explicit criteria are presented in Sec. Ill for the 
stability of each of the phases in the vicinity of bi-criticality. We construct the phase diagram 
at finite temperature and find excellent agreement with the experimental results of ReiclH, 
et. al. on the magnetic system LiHo x Yx-. x F&. The explicit solution for M > 2 is presented 
in the last section with a special emphasis on the Gabay-Toulouselll line. In section IV we 
analyze the possibility of replica symmetry breaking along the de Almeida- ThoulessS line 
in the M=l case. 



II. LANDAU ACTION 

Central to the construction of a Landau theory of the bi-critical region is the free-energy 
functional. For a quantum mechanical system, this is obtained by explicitly including in the 
partition function, 

(5) 

the time evolution according to a reference system. For this problem, the Hamiltonian for 
free quantum rotors 

describes our reference system. The perturbation 



Tr (e 



Tr 



-PHof 



exp 



H^dr 
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corresponds to the random magnetic interactions. The trace in the partition function is 
evaluated over the complete set of quantum rotor states. A primary hurdle in evaluating 
the partition function is the average over the random spin interactions. It is now standard 
to perform this averagei by replicating the spin system n times and using the identity, 

ln[Z] = lim ^ - 1 (8) 

to obtain the partition function, Z . We first must then evaluate (Z n ). Formally, the repli- 
cated partition function is defined through Eq. (§) by replacing H and Hi by their replicated 
equivalents 

n 

Hf = Y.H? (9) 

0=1 

where the superscript indexes the individual replicas and i = 0,1. Within the replica 
formalism, the average over the random interactions with the Gaussian distribution gives 
rise to a quartic spin interaction. This quartic spin interaction is easily reduced to a quadratic 
interaction by use of the Hubbard-Stratanovich transformation. Let us define the Fourier 
transforms 

S a (fc,r) = ^£S?(r) e ^ (10) 



of the local spin operator for each site and the corresponding transform of the nearest- 
neighbour interaction 

d 

J(k) = J^e k ' ri J = 2Jj2coskia (11) 

The replicated partition function now takes on the form, 

[Z n ] = Z£ J D^?DQe~ FeB [Q, (12) 
where the effective Free-energy 
F eS ^,Q] = Y / fdr^ l (k,T).(^(k,T)r+ £ Qf v (k,k',r,r')[Qf u (k,k',T,r')r 

a,k ^ a,b,k,k' 

(fexp \2 jyrJ2\fW%(k,r)S;(-k,r) 
+ E f f drdr'^2J(k)J(k')Q, u (k, k', r, r')S;(k, r)S b u {-k\ r') I ) 

a,b,k,k' J0 J0 I 



-In 



(13) 



is now a functional of the auxilliary fields Q and \1/ which appear upon use of the Hubbard- 
Stratanovich transformation to decouple the quartic spin term proportional to J 2 (k) and 
the quadratic term scaling as Jo(k), respectively. In writing this expression, we used the 
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Einstein convention where repeated spin (but not replica) indices are summed over, Zq = 
Tr[exp(— (3H )], and 



(14) 



The fields Q and \1/ play fundamentally different roles. The proportionality of the \l/ field 
to the mean of the distribution implies that this field determines the ordering transition 
(superconducting for M=2 or ferromagnetic for M=l). This can be seen immediately upon 
differentiating the free energy with respect to The self-consistent condition is that 



the order parameter for the ferromagnetic or superconducting phase within Landau theory. 
Likewise, differentiation of the free energy with respect to Q reveals that 



is the self-consistency condition for the Q matrices. For quantum spin glasses, it is the 
diagonal elements of the Q-matrix 



in the limit that |r — r'\ — > oo that serves as the effective Edwards- Anderson spin-glass order 
parameter!'! within the Landau theory. A disclaimer is appropriate here as D is non-zero 
in the spin-glass as well as the paramagnet phases. However, its behaviour is sufficiently 
different in the three phases: in the paramagnetic and ferromagnetic (superconducting) 
phases, as we will show, D{u) still has the form \J lo 1 + A 2 . The gap, A, vanishes at the 
transition to the spin-glass phase giving rise to the long-time behavior of D{j). 

Before we analyze the form of the gap, we must obtain the effective Landau action. The 
goal here is to obtain a polynomial functional of Q and \1/ from which a saddle point analysis 
can be performed. We proceed in the standard way and perform a cumulant expansion on 
the free energy. For the case of zero-mean, this procedure is documented quite closely in 
RSY. In addition to the terms containing the Q matrices studied by RSY, the non-zero mean 
case will contain powers of ty a as well as cross terms. The resultant action must contain, 
of course, only even powers of For an analysis of the bi-critical point, it is sufficient 
to retain quadratic and quartic terms in Terms of this kind are completely analogous 
to those derived previously by Doniachl. Of the cross terms, the simplest are of the form 
Wj^lQfy, ty^lQ™, and ^^QH. We have confirmed explicitly that retension of the 
latter two terms in which only one replica index occurs leads only to the renormalization 
of various coupling constants and minor modification of the phase boundaries near the bi- 
critical region. Hence, we do not consider such terms. Likewise, we do not retain terms of 
the form Q aa Q bb because this term vanishes in the limit n — > 0. 

We find then that the effective action 




Qf l/ (k,k\r,r') = (S;(k,r)S b u (k\r')) 



(16) 



(17) 
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dT l( iT2dT 3 ]T Qfu(x, n, T 2 )Ql C p (x, T 2 , T 3 )Q™{x, T 3 , 7i) 





d d x / dndr 2 £ T 2 )Q»t(^ ri, r 2 ) ■ ■ ■ 



a,6 



(18) 



contains three types of terms. 1) The terms that depend purely on the Q-matrices, as in 
the terms in the first curly bracket in Eq. (fHf ), have been studied previously by RSY in 
the context of quantum spin-glass/paramagnet transition. 2) The ^/-dependent terms in 
the second curly bracket give rise to the ferromagnet or superconducting phases. They are 
controlled by the parameter 7. 3) The competition between these two transitions is mediated 
by the last term in Eq. fll8[) . To re-iterate, other cross-terms exist-for example the term 
noted previously. However, such terms have no bearing on the critical region. It suffices 
then to truncate the action at the level of Eq. (0). The last term in Eq. fll8|) bears a strong 
resemblance to the term that appears in the Landau action for the zero-mean problem in 
the presence of a magnetic field, h. In fact, we can obtain the corresponding term from 
Eq. fll8|) by the transformation — > hS^iJ Kg jit. This mapping is not unexpected. In 
the non-zero mean problem, ordered spins create an effective "magnetic field" that acts as a 
source term for the Q-matrices. This analogy is particularly powerful and serves as a useful 
check as to the validity of our saddle-point equations. Further, this mapping is true even in 
the classical case. 

A word on the coupling constants is in order. The parametrization of the action in terms 
of the coupling constants k, t, and g was obtained by appropriately rescaling the fields Q 
and \I/ as well as the space and time coordinates. Fundamentally, g is a function of Ec and 
Jo, while k and t are functions of Ec and J only, u, v, and ( are related to a four-point spin 
correlation function evaluated in the long-time limit as discussed by RSY. In the zero-mean 
case, the phase diagram demarcating spin-glass and paramagnetic stability is determined by 
the parameter r which determines the strength of the quantum fluctuations. When these 
fluctuations exceed a critical value, that is r > r c , a transition to a paramagnetic phase 
occursi. In the problem at hand, two parameters, r and 7, determine the phase diagram. 
The coupling constant, 7, is directly related to the parameter Ec/ Jo and Ec/ J as well. The 
latter dependence arises as a result of the removal of the quadratic term, 




a,b 
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(19) 



by the transformations Q — > Q — C5 ab 5^{ri — r 2 ). From microscopic considerations, it 
follows that kH/2 = 1. This is an important simplification because we will show that the 
ferromagnet /paramagnet boundary is determined by the line 7 = 2A/k 2 £ = A. 



III. SADDLE POINT ANALYSIS: PHASE DIAGRAM 



u) = \ f dr^ a (k, r)e-^ T (20) 



In terms of the frequency-dependent order parameters 

-I 

(3 Jo 
and 

Qf u {k u k 2 ^L0 2 ) = j 2 j* dT X dr 2 Q a ^xM^T 2 y-^-^\ (21) 

the three phases in our problem form under the following conditions. For the paramagnetic 
phase near the bi-critical point, the saddle-point equations are satisfied when = but as 
RSY have shown 

Qf u {k,u; h uj 2 ) = (2n) d p6,J ab 5 d (k)5 UJl+UJ2fi D(u 1 ). (22) 

We have retained the ansatz for the Q-matrix used by RSY for the paramagnetic phase. 
Similarly, the ordering parameter also vanishes in the spin-glass phase, = and for the 
replica-symmetric solution, 



Qf u (k, u u u 2 ) = {2-K) d 5 d {k)5, v [pD(u l )5 ul+u>2 , 5 ab + /? 2 <Uo<Wg af> J . (23) 

The possibility of a replica broken-symmetry solution will also be explored by including 
terms of 0(Q A ) in the Landau action. In the ordered phase (that is ferromagnetic or super- 
conductor) both \1/ and Q are nonzero. As ferromagnet ism has generally been studied with 
a frequency-independent order parameter, we explore a static ansatz of the form 

u) = (27r) d 5 d (A0/?<L,o<WV> (24) 

for and 

fj, 



0%,(k, u 1} uj 2 ) = (2n) d 5 d (k)5 liu [PD(u l )6 IJl ^, 6 ab + (3 2 q5 UJlfi ^ 2 , S ab + /3 2 <U,o<W<7° b J (25) 

for the Q-matrices. Our explicit inclusion of q as the u — diagonal element of the Q- 
matrices implies that we can redefine D(u) = D(uj) + /3q5u,o- Consequently, we can set 
D(oj) = and assume that q ab is purely off-diagonal. The replica-symmetric solution corre- 
sponds to q ab = q for all a 7^ b. Initially, we will explore only this case. 
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A. M=l: Ferromagnetic order 




We now specialize to the M — 1 ferromagnetic crucial differences can occur with 

the analysis for M > 1. If we substitute Eqs. (|24D and (|25|) into the Landau action, we 
obtain a free-energy density of the form 

z ~ Trq 2 + Trf \ _ f3^_ f E ab q ab \ (2g) 
n n J Kgt \ n J 

Formally, the free-energy density as defined here is the disorder-averaged free energy per 
replica per spin component. The last term in Eq. ( f26"D arises from the cross term and 
generates off-diagonal components of the Q-matrix. For ip = 0, this term is absent and 
we generate the solution of RSY. We obtain the saddle point equations by differentiating 
Eq. ( p6|) with respect to q, q, D, and ip in the n — > limit. All of these quantities can be 
simplified once the appropriate gap parameter, 

A 2 = r + K(u + v) fg+^E^M] ( 27 ) 

is identified. The difference with the corresponding quantity in the work of RSY is the 
presence of the q term in the gap which includes the contribution arising from the cross 
term in the Landau action. From the derivative equation with respect to D(u), we find that 



D(u) = — Vuj 2 + A 2 . (28) 



K 



The constraint, A > is crucial for the stability of the free energy density and the minus 
sign in the above expression ensures that D(t) > 0. The saddle-point equations for q and q 
yield that 

(29) 



2ks-A 
and 

if) 2 A 

~ q = 2^R--p- (30) 

As expected, these equations form the basis for the bulk of our analysis and are identical to 
the RSY saddle-point equations with h replaced by ip. Despite this mapping, h and ip do 
serve different roles in the zero and non-zero mean theories. In the former, h is an external 
adjustable scalar quantity without any critical properties, whereas ip plays the role of an 
order parameter in the non-zero mean case and must be treated on equal footing with Q. 
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The saddle-point equation for if) 



if) |— 7 + A + (ip 2 \ = (31) 
implies that aside from the trivial solution if> = 0, the non-trivial solution corresponds to 

V> 2 = ^(7-A). (32) 

In obtaining this equation, we used the fact that K 2 t/2 = 1. If we specialize to low tem- 
peratures, that is low relative to the gap (T < A), the sum in the gap equation can be 
evaluated 

A 2 A 2 , A 2 



^EVOT = - + -m # (33) 

and the resultant self-consistent equation takes the form 

* 9 u + v . 9 , A 2 u + v ( 7 A ,„ , . 

A = r - r -— A hl # + ^(l- 1 )' (34) 

We have introduced 

r c =(« + ^ (35) 

and the cut-off A^ is determined by the energy scale at which zero-point quantum fluctu- 
ations become important. The natural cut-off for such fluctuations in this problem is Ec, 
the charging energy. 

The essential physics of the bi-critical point is contained in Eqs. ( |32"D and ([34]). From 
Eq. fl32|), it is clear that the ferromagnetic phase exists only when 7 > A. This ensures 
that if) 2 > in the ferromagnetic phase. The vanishing of the order parameter if) signifies 
a termination of the ferromagnetic phase. For 7 7^ and A ^ 0, the only line along which 
if> 2 = corresponds to 7 = A. If we substitute this condition into the gap equation, we find 
that the critical line separating the ferromagnet from the paramagnet at low temperatures 
(T <C A) is given by 

7 = ^= { % ) (36) 

\{u + v) ln^ 

and is depicted in Fig. (0). In obtaining this equation, we assumed that \r — r c \ 1 and 
7< 1. The bi-critical point corresponds to r = r c and 7 = 0. At this point, the gap vanishes 
as does if). The essential non-trivial nature of this result is the non-analytic dependence of 
7 on the quantum fluctuation parameter, r — r c . 

In the vicinity of the bi-critical point, distinct regimes partition the ferromagnetic phase 
that are determined by the magnitude of A and if). This behavior is shown in Fig. ([!]). The 
regions 01 and 02 are distinguished by their distance from the line 7 = A. The transition 
between 01 and 02 occurs when 
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Because r — r c is much less than unity, algebraic dependence in Eq. (p8[) of the form, 
(r — r c ) 3 / 2 , implies that the 02 region is quite narrow. Within this narrow region, the gap 
is well approximated by Eq. (|36D and the value of ip is given by 

^ 2 = 7h- L 47r(r ~ r i ) | region 02 (38) 
[u + v) In 

In the bulk of the ferromagnetic phase, region 01, the 7-dependent terms must be retained 
in Eq. (|34D to accurately describe the gap, A = 7 — (ip 2 , with the ordering parameter given 
by 




,2 2 9 
iP = — — 

U + V 



— — 7 In l(r-r c ) 

Z7T 7 



region 01 (39) 



Hence, at the point r = r c , we find a strong non-analytic dependence, ip oc 7 3/ ' 2 y / ln A w / 7 , 
on the coupling constant 7. This is particularly important because it signifies that even at 
the mean-field level, deviations from the standard square-root dependence are present in the 
current theory. This result is fundamentally tied to the logarithmic dependence induced by 
the frequency summations and is caused by the zero-point quantum fluctuations. However, if 
we extrapolate our results to the regime, r—r c « 0(1), the second term in Eq.(p9|) dominates 
and we do recover that ip oc 7 1 / 2 in agreement with the expectation from standard mean-field 
Landau theories. Paramagnetic and spin-glass behaviour obtain whenever 7 < A. In this 
regime, the non-trivial solution for ip no longer holds and ip — is the only valid solution. In 
this limit, our solution for A is identical to that of RSY and all of their results are recovered. 
For example, consider the free energy density (in units of K 2 t/2 = 1), 

£ = (- + u + v \ _ A K r - r e) _ ( r - r cf l 2 (r - r c ) _ 7* ^ _ _ _ 

n 2tt V6 8tt / 4tt 4(u + v) 2(u + v) 8tt ^ 7 

8tt 7 2(u + v) K J 



in region 01. This quantity is obtainable from Eq. ( p6[) once the saddle point solutions for 
region 01 are used and only the leading terms in 7 and the cross-term are retained. When 
7 = 0, this expression is identical to that of RSY in the spin-glass phase. As anticipated, the 
leading 7 dependence in the free-energy density is non-analytic. This behaviour originates 
from the non-analytic behaviour of the order parameter ip on 7. We will see that this 
non-analyticity does not survive for M > 1 above and below the Gabay- Toulouse line. 

Extending these results to finite temperatures, T ^> A, simply requires the evaluation of 
the frequency summation over frequencies in the gap equation for T ^> A. Using the result 
in Eq. (2.13) of RSY, we obtain a self-consistent condition for the gap 

A 2 = r - r c (T) - (u + v)TA - ^A 2 In ^ + - A) (41) 
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in the regime T ^> A where we have defined r c (T) = r c — (u + v)ttT 2 /3. Recall that the 
transition between the ferromagnetic and paramagnetic phases occurs when 7 = A. Hence, 



in the units chosen here, this condition simplifies to 7 = A(T). For T <C yr — r c (T), the 
boundary for the paramagnetic/ferromagnetic state remains unchanged from the T = 
results discussed above. However, for T ^> A, two distinct regimes 



r , rc( - T ^ \ r — r c (T) <C T region 03 

7 = A(T)H ^ , ' n/ , (42) 

"" ' V r - r c < T region 04 



31n%^ 



emerge depending on the magnitude of the thermal fluctuations. These regimes are depicted 
in Fig. (0a). The crossover between these two regions occurs when 7 oc y/r c — r. In Fig. (0a), 
r c — r > 0. The temperature 



Tn 



V ^4 («) 

\ ir(u + v) V 7 



is denoted explicitly in Fig. (]|a) as this is the lowest temperature at which region 03 
obtains. Immediately below region 03 where 7 — A (T) oc (r — r c (T)) 2 /T and to the right 



of 04 where 7 — A (T) oc T 3 / yln A w /T, a transition to a new region occurs in which the 
gap takes the form, 



7 



u + v 



— In — + (u + u)7 T - (r - r c (T)) 7 



7 - ^ 2 {T) (44) 



In this region, denoted as 05 in Fig. (0a), as well as in regions 03 and 04, classical thermal 
fluctuations dominate. These regions can be construed as being quantum critical. Further 
away in region 01, the ferromagnetic phase is impervious to thermal fluctuations. The 
crossover to this regime occurs when 7 w 0(T). This partition is the dashed line separating 
region 05 from 01. In Fig. (@b), the corresponding phases are shown for r — r c > 0. The 
key difference with the r — r c < regime is the absence of the spin-glass phase. 

Experimentally, the phase diagram has been measured for the random Ising spin system 
LiHoxYi^xF^ at finite temperature. This system possesses all three phases discussed here. 
It then serves as a bench mark test of the phenomenological theory we have developed. 
While the overall features of the experimentally-determined phase diagram are similar to 
that shown in Fig. (0a), it is worth looking closely at the form of the boundaries between 
the three phases. Particularly striking in the experimentally-determined phase diagram^ 
is the close to linear dependence of the PM/FM phase boundary away from the bi-critical 
region but a non-linear dependence on the doping level in the vicinity of the bi-critical 
region. This dependence mirrors closely the behavior of the PM/FM finite temperature 
phase boundary shown in Fig. (0a). While a quantitative comparison cannot be made 
because of the phenomenological nature of the coupling constants used in this model, the 
agreement with experiment is sufficiently striking and serves to justify the applicability of 
the model used here. 
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B. M > 1 



We consider now explicitly M > 1. For the problem at hand, the ordered phase for 
M = 2 corresponds to a superconductor. Analogous isotropic solutions can be obtained for 
M > 1 with the transformations u+v — > u+Mu. However, because non-zero mean generates 
spontaneously an effective magnetization, there exists a possibility that the different spin 
components of the replica Q-matrices might acquire fundamentally different values as first 
proposed by Gabay and Toulouse^. In the zero-mean case, this happens only when a 
magnetic field is present. However, in this case, the Gabay- Toulouse (GT) line exists for all 
M > 1 as a result of the spontaneously-generated magnetization. 

To explore the possibility of a GT line, we must generalize the ansatz for the Q-matrices 
to explicitly break the symmetry between the spin-components of Q. The simplest way of 
doing this is to divide the spin components of the Q-matrix into longitudinal, \l — v — 1 
and transverse, \i — v ^ 1 sectors. Hence, in Eq. (P5|), we introduce the parameters, 
qL, and D l {uj) for the longitudinal fj, — v — 1 component and qft, Qt, and D T {ui) for the 
transverse components fi > 1. At the replica-symmetric level, both g£ b and are constants 
independent of the matrix label, ab. We will call these constants, qi and qr, respectively. 
The resultant expression for the free-energy 



n (5Kt 



(3nt 



+ (m - 1)^ - ^ e (^iH + (m- m 



u 

+ 2t 



rqx k 
~Kt 3/3t 

2 



+ (M — 1) 



<?£ + 75 E D t( u ) 

P LU^O 



.Of. 



(ql ~ Zhql + 2g£) — (M — l)^(g| - 3g T <# + 2 4 



3t 



{qL - qh) + — -7^ + "TT 



1 

2# 



(45) 



is a generalization of Eq. (|26| ) to an anisotropic system. The explicit factor of M — 1 arises 
from the separation into transverse and longitudinal components. 

If we approach the GT line from below, we find that the relevant saddle-point equations 

are 



A£ = r + [u + v) K [q L + - E ^Ml + ( M - !)««[«r + ^ E ^M] 



K 
1, 

-—a; 
n 



(46) 
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Ql 

Ql 




2ngA L 

r al 

2ngAi k/3 



q T = q T 



E 



CO 



u + (M-l)v \/3 



E 



A| 




which are a direct generalization of the M = 1 equations to the anisotropic system. We are 
particularly interested in the solution in the 7-r plane where qx = 0. This demarcates the 
GT line. Below the GT line, q~T = Qt 7^ 0, while above, the transverse replica off-diagonal 
component of Q vanishes. Note this state of affairs does not occur unless ip ^ 0. If we 
substitute the non-trivial solution for ip into the expression for A|, we find that within 
logarithmic accuracy at zero temperature, we recover the result obtained previously for 
M = 1 but with u + v — > u + Mv. The phase diagram hence is identical to that shown in 
Fig. ([!]). However, a new region, 01, appears. This is illustrated in Fig. (0a). To find the 
line demarcating this region we must solve for the transverse replica off-diagonal component 
of Q. After several manipulations of the set of equations in Eq. (0), we find that 



1 



q T 



u 



k{u + Mv) 

If we use the fact that A^ w 7, we find that the GT line occurs when 



7 



u 



(47) 



(48) 



The phase diagram depicting this line at zero temperature is shown in Fig. (0a). In the 
region labeled 01, qr 7^ and qi 7^ 0, whereas in 01 only q^ 7^ 0. Hence, we have identified 
the zero-temperature GT line. At finite temperature, the generalization of the Eq.(|48|) is 
simply 



7 



\ 



r c — r — (u + Mv) 



ttT 2 



(49) 



Hence, the GT line is now a surface in the 7, r, and T space, a slice of which is shown in 
Fig. (|b). At the point 7 = 0, we recover the isotropic result that 



Qt = Ql 



r c (T) 



k(u + Mv) 



(50) 



in the spin-glass phase. 

We now consider the region above the GT line. This region was not analyzed by RSY. 
However, this region is of considerable interest because although q? vanishes in this region, 
the non-zero transverse component of the order parameter becomes gapped, significantly 
different in character from the longitudinal one. This can be seen immediately by setting 
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Qt = Qt = in Eq. fl45|) and differentiating with respect to Dt{uj). From this operation, we 
find that contrary to the ungapped transverse component of Q below the GT line, 

D T (co) = --(lo 2 + A T )^ (51) 



K 




with 



above the GT line. The corresponding expression for At is easily obtained from the first 
equation in Eq. ( pEB]) by setting = 0. The ferromagnetic phase above the GT line 
can be divided into two regions, 01 and 02, which are now different with respect to the 
relative magnitudes of At and At- In the region 02 which is completely analogous to the 
corresponding region in the M — 1 case, At and At are almost equal and are given by 
Eq. ( [36]) with u + v — > u + Mv. The condition for crossover to 01 in which At and At 
are somewhat different in magnitude is given by Eq. (|37[) . This conclusion is reached by 
manipulating the system of equations in Eq. ( J4"6j ) with q T = <?t — and D t {uj) given by 
Eq. (|51"D . The resultant expression, 

vA\ - (u + v)A T = -ur + < U + Mv ^ ^ ^2 + ^ (53) 

P OJ 

which is valid at any temperature, contains both the transverse and longitudinal gaps. 
Within the approximation that At pa 7, we obtain that at zero temperature in 01, 



47r(t>7 + u(r — r c )) 
u(u + Mv)\n- A ~ 



(54) 



(vy 2 +u(r— r c )) 

From this equation it is immediately clear that At is logarithmically smaller than At- Also, 



we easily recover the GT line, 7 = yu(r c — r)/v , simply by solving At = 0. 

At finite temperature, we can formally distinguish three limiting cases: 1) At 3> T and 
At > T, 2) A L > T and A T < T, and 3) A L < T and A T < T. The regime A L < T 
and At 3> T does not exist as At is always less than At- Case (1) is identical to the 
T = limit whereas cases (2) and (3) are the high-temperature limit with respect to At- 



To describe (2) and (3), we calculate the sum in Eq. (|53D in the high-temperature limit with 



the approximation that At ~ 7. Within logarithmic accuracy, the resultant equation 

iH 2 - u(r c (T) -r)= u(u + Mv) \ n ^ + T A T ^j (55) 

is similar in structure to Eq. (JO]). Consequently, within cases (2) and (3), two distinct 
regimes denoted by 05 and 05' and 06 and 06', respectively, arise. In the regions super- 
scripted with a prime, two conditions hold 

- (r c (T) - r) < T, A T = — - u{r ^ - r) 05 and 06 (56) 
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Contrastly, in the unprimed regions, 



v 7 2 — (r c - r) < T, A 3 



2ttT 2 



31n^ 



05' and 06'. 



The 



transition 



between 



the 



primed 



and 



(57) 



unprimed 



regions occurs when ^Jvy 2 — u{r c (T) — r) w T/lnA^/T, whereas the transition between 
06' and 01 occurs when 



't>7 2 — u(r c — r) fa T. The difference between regions 05 and 
05' and regions 06 and 06' is that T <C Ax in the latter whereas the opposite is true in 
the former. We have taken particular care in distinguishing the primed from the unprimed 
regions because they imply that the GT transition is identical in structure to the ordinary 
paramagnet/spin-glass transition described by RSY. Except, only the transverse component 
of Q is affected at the GT transition. In fact, regions 05' and 06' are quantum critical with 
respect to the GT transition while the temperature dependence of the transverse component 
is 'classical' in regions 05 and 06. In 01, thermal fluctuations are subservient to quantum 
fluctuations for both transverse and longitudinal components of Q. 

In each of these regions, a key question that can be addressed is how does the transverse 
gap renormalize ip and Ax,. Consider first the regime 7 ^> T. In this regime, we find that 



A/ 



7 - C^ 2 (T) = 7 



u + Mv 



2/T 7 



(M-l> 7 2 
T)1 u{u + Mv) T 



(58) 



It follows immediately from Eq. (|54|) that the transverse contribution to the longitudinal 
gap is logarithmically small. As a result, the order parameter ip is well described by its value 
given by Eq. ( |39|) with u + v — > u + Mv. Consider now the high-temperature regime. In this 
case, Ax is exactly given by Eq. (pH|) plus the term proportional to A 2 in Eq. Q5B|). However, 
we checked explicitly that in every sub-regime, the correction due to Al| is sub-dominant to 
the leading terms. Consequently, the order parameter, ip, and Ax are both unrenormalized 
by the transverse gap to the leading order. 

Using the solutions delineated in Eq. ((46[ ), we can calculate the free-energy density below 
and above the GT line. Above and below the GT line, we find that the free-energy density 
at zero-temperature up to the leading 7— dependent terms, 



n 



Fq{u + v 



u 



Mv,C) 



(M - l)vA l^ 
Attu 



+ 



7 2 ( r 



2u(u + Mv) 



+ 



u 



(59) 



contains the standard analogous contribution from the M = 1 analysis as well as qx and 
Ay-dependent terms arising from Eq. fl45|) . It is the contribution from the latter terms that 
results in a suppression of the 7 4 lnA a; /7 as the leading 7— dependent terms when r = r c . 
At the bi-critical point, we find that in contrast to the M=l case, the leading term in the 
free energy density is analytic in the coupling constant, 7. This term is of the form, A 2 7 2 . 



IV. REPLICA-SYMMETRY BREAKING 

The requisitelll for a replica asymmetric solution within Landau theory of spin glasses is 
the presence of a Q Ath term in the Landau action. We specialize to M — 1 for simplicity. To 
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facilitate such an analysis, we must extend the cumulant expansion of Eq. ([13]) to the next 
order in perturbation theory. As RSY have performed such an analysis for the spin-glass 
phase, we focus on the ferromagnetic case. While several types of fourth-order terms occur, 
the most relevant is of the form, 



yi 

Gt 



J d d x J dndT 2 £ [Q ab (x, n, r 2 )] 4 (60) 

ab 



This term will give rise to a (g afe ) 4 contribution to the free-energy. Our focus is the resultant 
change in the free energy 



^ = -^(7V + 8,ZV)-g S9 M-^ § (^. (61) 



The presence of the ip 2 term suggests that within the space of ultrametric functions^ 
q(x) on the interval < x < 1, we should choose an ansatz for q, 

!q < s < s 
If so < a < 8 X (62) 
qx si < s < 1 

that has two distinct plateaus. This insight is based on an analogy with the replica broken- 
symmetry solution in the presence of a magnetic field. From continuity, we must have 
that go = ^/9so/2yi and q± = K(3si/2yi. The constants so and s\ can be determined from 
the saddle-point equations for AjF. Upon differentiating with respect to qi, we find that 
q = q\ — y\q\j ' k(3. The corresponding equation for q provides a relationship 

( 3 ^ \ k 

between q and ip. Replica symmetry breaking occurs when q Q < q\. To leading order in 
temperature, we can approximate q\ ~ q, where q is given by 

Q = — — r- (64) 
Hence, the boundary demarcating the replica-symmetric solution is determined by 



4 



= 3 A (65) 

which we obtain upon equating q and q\. If we use Eq. (PD|) for ip which is valid in region 
01 and use the fact that in this region A pa 7, we obtain 

2?/i(r c — r) 2 

1= ^g(u + vf r<rc (66) 

as the condition for replica symmetry breaking at low temperatures, T 7. This condition 
for replica-symmetry breaking extends continuously to the spin-glass phase agreeing with the 
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work of RSY. The phase diagram illustrating the replica-broken symmetry region is depicted 
in Fig. (|^). At finite temperature, we use Eq. ([Hj) for ip and obtain a generalization, 



7 



2 2/i I" rp i r c -r] 2 
7T + 



3n 2 g 



u + v 



r < r c , (67) 



of the replica-breaking condition which extends smoothly over to the zero-temperature con- 
dition. 

To estimate the strength of the replica-symmetry breaking, we define the effective broken 
ergodicity order parameter, 

A, = gi- I q{s)ds. (68) 
If we substitute the expression for q(s) and integrate, we obtain 

\ = —(Ql(T)-q 2 (T)). (69) 



K 



Because both q and q\ are finite at low temperatures, A q — > as T — > 0. The weakness of the 
replica-symmetry breaking in the ferromagnetic phase is in accord with the weak-symmetry 
breaking found by RSY in the spin-glass phase. Implicit in the replica-broken solution 
in Eq. (|6"2]) is the presence of many degenerate energy minima in the energy landscape. 
The weakness of replica symmetry-breaking at low temperatures in this model within the 
ferromagnetic phase suggests that the ferromagnetic phase is energetically homogeneous. 



V. SUMMARY 

We have constructed here a Landau theory near the bi-critical point for a ferromagnet, 
spin-glass, and paramagnet. The analogous analysis was also performed for M > 1 in 
which the ordered phase for M = 2 is a superconductor. All transitions were found to be 
second order in contrast to the assertion by Hartman and Weichmani who claimed that 
the FM/PM transition was first order in the spherical limit. A key point of our analysis is 
the formal equivalence between the role of a non-zero mean and the presence of a magnetic 
field in the zero mean problem. This observation is equally valid for classical systems. 
Resilience of the ferromagnet against thermal fluctuations occurs when T < 7, where 7 
is the coupling constant that ultimately determines the rigidity of the ferromagnet phase. 
Additional features of our analysis that are particularly striking are 1) the non-analytic 
dependence of ip on 7, namely ip oc 7 3//2 / In A^/7, near the bi-critical region, in the vicinity of 
the bi-critical point and 2) the subsequent leading non-analytic dependence of the free-energy 
density in region 01 on 7 for M = 1. We showed that this behavior ceases for M > 1 below 
and above the GT line. The reason underlying this difference with the M = 1 case is the 
presence of M — 1 components of Q that yield leading analytic contributions of order 0(7 2 ) 
to the free energy density. An additional feature which our analysis brings out for M > 1 
is the similarity between the Gabay- Toulouse transition with the spin-glass/paramagnet 
transition with zero mean. This similarity has been noticed previously in the context of 
classical spin-glasses0. A key surprise found in the analysis of the GT transition is the sub- 
leading depdence of ip and on the transverse gap and q?. This suggests that the GT line 
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should have only weak experimentally-detectable features in the superconducting phase, for 
example. The excellent agreement observed with the experimental results on LiHo x Yi^ x F^ 
for the case of M = 1 is encouraging that similar agreement will be found with experiments 
on the analogous superconducting systems. 

While we referred to the point of intersection between all the phases as being bi-critical, 
it is in fact multi-critical. This state of affairs obtains as a result of the presence of replica 
symmetry breaking and the Gabay- Toulouse instability. As in the M — 1 spin-glass case, 
we also showed that the non-ergodicity parameter is linear in temperature illustrating the 
weakness of replica symmetry breaking in the ferromagnetic phase. Though we did not 
treat explicitly replica-symmetry breaking for M > 1, such symmetry breaking is expected 
in this case as well when quartic terms are included in the action. In a future study, we will 
extend this analysis beyond mean field and report on the renormalization group analysis of 
the 'bi-critical' region. 
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FIGURES 



FIG. 1. Zero-temperature phase diagram demarcating the regions of ferromagnetism (FM), 
paramagnetism and spin-glass behaviour for M=l. The parameters 7 and r are determined by the 
dynamical quantum fluctuations and the static disorder. The curve separating the ferromagnet 
from the paramagnet scales roughly as 7 ~ y/r — r c up to logarithmic factors. Regions 02 and Ol 
are distinguished by the magnitude of the spin- wave gap, A, as well as the magnitude of the order 
parameter ip, which is non-zero only in the ferromagnetic phase. The transition between these two 
regions occurs when 7 — 7 C ~ (r — r c ) 3//2 , where 7 C and r c are determined by Eqs. ( |36| ) and (|35|), 
respectively. 

FIG. 2. M=l finite-temperature phase diagrams for a) r c — r > and b) r — r c > 0. Regions 
03 and 04 lie close to the paramagnet (PM) boundary and hence have 7 ~ A(T). Thermal 
fluctuations dominate in region 05 as well as in 03 and 04. In region Ol, thermal fluctuations are 
negligible. The key difference between noted when r c — r changes sign from positive to negative is 
the absence of the spin-glass (SG) phase for r > r c . 

FIG. 3. a) Zero and b) finite-temperature phase diagrams for M > 1. The criteria for distin- 
guishing regions Ol and 02 are identical to the M = 1 case but except u + v — ► u + Mv. The 
GT line separates regions Ol from Ol. This line terminates at 71 = \/u(r c — r) /v. Below this 
line both the transverse and longitudinal components of the replica off-diagonal components of the 
Q-matrix are non-zero. The location of the GT line is given by Eq. (f48|). b) Finite temperature 
GT line as determined by Eq. (|49|). 

FIG. 4. Phase diagram illustrating the distinct regimes that occur at finite temperature above 
the GT line denoted with solid circles. Regions 03 and 04 are as described in Fig. (|2|a) previously. 
Regions 05' and 06' quantum critical with respect to the GT transition. The difference between 
regions 05 and Oh' and regions 06 and 06' is that T <C A^ in the latter whereas the opposite is 
true in the former. 

FIG. 5. a) Phase diagram illustrating replica-symmetry breaking for M=l in the 7 — r 
plane. The broken symmetry region is denoted RSB. Regions Ol and 02 are as before, b) 
Replica-symmetry breaking region (shaded region) at finite temperature. The finite temperature 
criterion is given by Eq. ([37]). The point 72 = 2y\(r c — r) 2 / \3k 2 g(u + v)) is determined by the 
replica-symmetry breaking condition specified in Eq. (|67]). 
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